Inverse problem of variational calculus for nonlinear evolution equations 
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We couple a nonlinear evolution equation with an associated one and derive the action principle. 
This allows us to write the Lagrangian density of the system in terms of the original field variables 
rather than Casimir potentials. We find that the corresponding Hamiltonian density provides a 
natural basis to recast the pair of equations in the canonical form. Amongst the case studies 
presented the KdV and modified KdV pairs exhibit bi-Hamiltonian structure and allow one to 
realize the associated fields in physical terms 
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^ 1. Introduction 
a- 

' In the calculus of variations one is concerned with two types of problems, namely, the direct and inverse problems. 

y—i , The direct problem is essentially the conventional one in which one first assigns a Lagrangian and then computes 

^-H ■ the equation of motion through Euler-Lagrange equations. As opposed to this, the inverse problem begins with the 

[ equation of motion and then constructs a Lagrangian consistent with the variational principle The object of the 

HH ' present work is to derive an uncomplicated method for the Lagrangian representation of nonlinear evolution equations. 

C/^ ' We shall see that our results for the Lagrangian densities provide a natural basis to recast these equations in the 

. Hamiltonian form 0, Q . 

• ^ ' Studies in the Hamiltonian structure of nonlinear evolution equations are based on a mathematical formulation that 

, does not make explicit reference to Lagrangians Q • We feel that the Lagrangian formulation of these equations should 
be quite interesting because Lagrangian densities, via the Legendre map, will give us a direct route to construct the 

, expressions for Hamiltonian densities that characterize the equation of Zakharov, Faddeev 0] and Gardner [3]. To 

^ ' gain some weightage for the physical and mathematical motivation of our work we proceed by noting the following. 

I Let P[v] — P{x, u*^"^) G be an r-tuple of differentiable function. The Frechet derivative of P is the differential 

■ operator Dp : A"^ ^ defined by 

o 

§: DpiQ)^^\e=oP[v + eQ[v]] (1) 



o 



for any Q G A'^. If P = ^ P[u]Dj, Pj G ^ is a differential operator, its adjoint D* is given by 
J 



D^^Y.^-D)j.Pj. (2) 



J 



I: 

■ a Lagrangian representation only if Dp is self-adjoint. When the self-adjointness is guaranteed, a Lagrangian density 
C for P[u\ can be explicitly constructed using the homotopy formula C[v\ = vP[\v\d\. A single evolution equation 
is never an Euler-Lagrange expression. One common trick to put a single evolution equation into the variational form 
is to replace u by a potential function w with u = —Wx- This yields Wxt = P[wa;].The function w is often called the 
Casimir potential. For many nonlinear evolution equations the Frechet derivative of P[wx] is self adjoint, while that 
of P[u\ is not. In view of this, the Lagrangian densities for single evolution equations are always written in terms of 
partial derivatives of w. There are some equations for which even P[wx] is not self-adjoint. These are often referred 
to as non-Lagrangian. The well-known Burgers equations and nonlinear evolution equations with nonlinear dispersive 
terms [6j] serve as typical examples of non-Lagrangian equations. 

Keeping the above in view we shall work out a Lagrangian representation without taking recourse to the use of 
Casimir potentials. In section 2 we introduce a suitable associated equation for an additional field variable v{x, t) and 
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use it to write the action principle for any nonlinear evolution equation. We express the Lagrangian density in terms 
of u{x,t), v(x,t) and their partial derivatives. We show in section 3 that the coupled set of equations as introduced by 
us form a Hamiltonian system. We also cite examples in which the associated fields admit simple physical realization. 
In section 4 we summarize our outlook on the present work and make some concluding remarks. 



2. Lagrangian representation 

The basic philosophy we use here to construct expressions of Lagrangian densities of nonlinear evolution equations 
has a rather old root in the classical- mechanics literature. For example, as early as 1931, Bateman |71] allowed for an 
additional degree of freedom to bring the equation of motion for the damped harmonic oscillator within the framework 
of action principle. The Bateman Lagrangian 

L ^ xij +^{xy - yx) ~ uj'^xy (3) 

for the one-dimensional linearly damped harmonic oscillator 

x + ^x + Lo^x — Q^ X — x{t) (4) 

has a 'mirror image' equation 

y~'iy + w^y = 0, y = y{t) (5) 

for the associated coordinate y{t). Here 7 is the co-efficient of friction and cj, the natural frequency of the oscillator. 
The overdots stand for derivatives with respect to time t. Understandably, the complementary equation in (5) 
represents a physical system which absorbs energy dissipated in the first. Interestingly, Bateman |7|] regarded a 
dissipative system as physically incomplete such that one needs to bring in an additional equation to derive the 
original one from an action principle. Thus it remains a real curiosity to envisage a similar study in the context of 
classical fields and look for the Lagrangian representation of nonlinear evolution equations. 

Any nonlinear evolution equation that has at least one conserved density p[u\ can be written in the form 

d 

-I- — = 0, u^u(x,t). (6) 

This single evolution equation is non-Lagrangian. When written in terms of the Casimir potential, the equation 
resulting from (6) may be either Lagrangian or non-Lagrangian. However, we can make use of an elementary lemma 
to get a Lagrangian representation of (6). 

Lemma 1. There exists a prolongation of (6) into another equation 

vt + ^{p[u]v^) ^ v^v{x,t) (7) 
ou 

with the variational derivative 

— = y(-i)^—-^ uk ^— (8) 

5u ^ dx^ dukx ' ^ dx'^ 

k—O 

such that the system of equations follows from the action principle 

6 [ £ dxdt = 0. (9) 



Here C stands for the Lagrangian density. 

Proof. For a direct proof of the lemma let us introduce C in the form 



C =^ ^{vut - uvt) - p[u]vx. (10) 



From (9) and (10) we obtain the Euler-Lagrange equations 



dt ^ dvt Sv ^ ^ 
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and 



dt dut 



SC 

Su 



(12) 



Using (10) in (11) and (12) we obtain (6) and (7) respectively. 



We shall now consider two examples of physical interest and apply the rule in (7) to construct the associated 
equations. We first focus our attention on the Korteweg-de Vries (KdV) equation 



Ut + 6uUx + MSx = 0. 



(13) 



The KdV equation represents the prototypical nonlinear evolution equation that was first solved by the inverse spectral 
transform method 8]. From (6) and (13) we see that for this equation 



p[u] = 3?i^ + U2x- 

Using (14) in (7) we get the associated equation 

The corresponding Lagrangian density as obtained from (10) reads 

£ = ^{vUt - UVt) - (3u^ + U2x)Vx- 



(14) 
(15) 

(16) 



It is easy to verify that (16), when substituted in (11), reproduces the KdV equation. The second example of our 
interest is the so-called modified KdV (mKdV) equation given by 



Ut + 6u Ux + u^x = 0. 



(17) 



Equation (13) and (17) are connected by Miura transform and as with (13), (17) can also be solved by the inverse 
spectral method 0]. The mKdV equation appears in a number of applicative contexts including description of Alfven 
waves in a collisionless plasma. The associated equation for (17) is obtained in the form 



Vt Vx+ V^x = 

with the Lagrangian density given again by (10). For the niKdV equation 



p[u] = 2u^ 



U2x- 



(18) 
(19) 



Results similar to those for KdV and mKdV equations can also be written for other nonlinear evolution equations 
which can be expressed in the form (6). We give in Table I, the results for a number of evolution equations which are 
often believed to be non-Lagrangian. 

Table I. Associated equations for a few important nonlinear evolution equations 



Evolution equation 


Conserved density 


Associated equation 


Burgers2: 






Ut - U2x - 2uUx = 


-{Ux + U^) 


Vt + V2x - 2uvx = 


Burgers3: 






Ut - U3x - Su'^Ux- 


-{U2x + + 3uMx) 


Vt + 'iuV2x — 'iu^Vx 


■iuU2x — 3u^ = 




-V3x = 


KdV-Burgers: 






Ut + UUx - VU2X + 


\u^ - VUx + flU2x 


Vt + UVx + VV2X + 


^J.U3x = 




AiW3a; = 


FNE3: 






Ut + Su'^Ux + 6UxU2x + 


U^ + 2ul + 2uU2x 


Vt + iu^Vx + 2uV3x = 


2uU3x = 






FNE5: 






Ut + f3i{u^)x + P2iu^)3x 


PiU^ + I32{u^)2x + P3{u'^)ix 


Vt + 2PiuVx + 2I32UV3^ 


+!33{u^)5x = 




+2/33UV5X = 
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The first two equations in the Table are due to Burgers. These are dissipative and do not support sohton solutions. 
However, both of them are useful in the study of acoustics and shock waves [l^I- In the recent past two of us 
studied the equations in the Burgers hierarchy and sought a Lagrangian representation in which the appropriate 
equations were expressed in terms of Casimir potentials. However, the results represented here have the obvious 
virtue of simplicity and directness. The compound KdV-Burgers equation describes wave propagation in which 
the effects of nonlinearity, dissipation and dispersion are all present. We believe that the Lagrangian representation of 
these equations using associated equations is quite interesting. The KdV and mKdV equations are quasi-linear in the 
sense that the dispersive behaviour of the solution of each equation is governed by a linear term giving the order of 
the equation. In contrast to this, the third and fifth order equations (FNE3) and (FNE5) ^ in the table are nonlinear 
partial differential equations with nonlinear dispersive terms. These are, therefore, fully nonlinear evolution (FNE) 
equations. The solitary wave solutions of these equations have a compact support. To the best of our knowledge 
no Lagrangian representations of FNE3 and FNE5 have yet been found. There have, however, been attempts [l^ to 
introduce Lagrangian system of FNE equations which support compacton solutions. 

3. Canonical structure 

Zakharov and Faddeev 0] developed the Hamiltonian approach to integrability of nonlinear evolution equations in 
one spatial and one temporal (1 + 1) dimensions and Gardner ^21, in particular, interpreted the KdV equation as a 
completely integrable Hamiltonian system with dx as the relevant Hamiltonian operator. In this context we introduce 
the following lemma for the Hamiltonian structure of the coupled equations introduced by us. 

Lemma 2. The Hamiltonian density Ti constructed from (10), can be used to express (6) and (7) in the Hamiltonian 
form 

Vt = J?^ (20) 



with 



and the symplectic matrix 



Sr] 



V-il) (21a) 



Proof. Using the Legendre map the Hamiltonian density for the Lagrangian in (10) is obtained as 

n = p[u]vx. (22) 
From (20), (21) and (22), equations in (6) and (7) follow immediately 

A significant development in the Hamiltonian theory is due to Magri |14i who realized that completely integrable 
Hamiltonian system have an additional structure. They are bi-Hamiltonian systems, i.e., they are Hamiltonian with 
respect to two different compatible Hamiltonian operators. We have found that both KdV and mKdV can be recast 
in the bi-Hamiltonian form 

■qt = Ji^— ^ -h^—- (23) 
6rj or] 

The appropriate results for Hamiltonian operators and Hamiltonian densities are given by 

jKav/mKav ^ (24) 

jKdv _( {dl + 2udx + 2dxu)d-^ \ 
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-y-idl + Au") 2uv,d-^ )' ^^^^ 

Hi = uv^, (27) 

- (3^2 + U2,)v, (28) 

and 

^™w^(2u3_^U2.)«x. (29) 

In solving the inverse variational problem for nonlinear evolution equations we coupled the field variable u{x,t) of 
some given equation with the field variable v{x,t) of an associated equation such that the system follows from the 
action principle with a prescribed form of the Lagrangian density as given in (10). Admittedly, one of our tasks in this 
work will be to study the nature of v{x,t) when the original field variable u{x,t) admits simple physical realization. 
Keeping this in view, we focus our attention on the KdV and mKdV equations which support soliton solutions. It is 
well-known that the solutions of (13) and (17) are given in the general form 

u^^'^ix, t) = A{k)sech^{kx ~ 4kh) (30) 

^rnKdVf^^^ t) = B{k)sech{kx ~ kH) , (31) 

where k is the wave number for a single bound-state energy of the potential that characterizes the spectral problem 
in solving the evolution equation [8]. Here A{k) and B{k) represent the amplitudes of the bright solitons in (30) and 
(31). Understandably, both solutions u^'^^{x,t) and u^^'^^ {x,t) are centered at x = 0. The KdV solution moves to 
the right with speed 4fc^ while the mKdV solution moves in the same direction with speed fc^ only. One can verify 
that the associated fields corresponding u^''-^ {x,t) and u"^^'^^ {x,t) are given by the dark soliton solutions 

^Kdv/mKdv^^^ t) = tanh(fca; - 4kh). (32) 

From (30) and (32) it is clear that both soliton solutions of the KdV pair move with equal speed. As opposed to this, 
comparison of (31) and (32) reveals that the speed of the dark soliton solution for the mKdV pair is four times the 
speed of the bright soliton solution. This is not immediately clear to us and deserves extensive numerical study for 
further clarification. 

4. Concluding remarks 

Nonlinear evolution equations are not directly amenable to Lagrangian representation. We have established that if 
(6) does not admit a direct analytic or Lagrangian representation, then there exists an auxiliary or associated field 
which helps us treat the original evolution equation within the framework of the action principle. The method is quite 
general and works for both integrable and nonintegrable equations. We have dealt with equations in which the effects 
of nonlinearity, dissipation and dispersion are all present. 

It is a well-known belief that there is no exact method for applying variational principle to dissipative systems. In 
view of this, studies in Lagrangian and Hamiltonian mechanics of nonconservative systems are still regarded as an 
interesting curiosity. In the context of point mechanics Riewe [15i] used the method of fractional calculus to write a 
Lagrangian for (4) without taking help of the additional equation in (5) . Kaup and Malomed [l3| sought an application 
of the variational principle to nonlinear field equations involving dissipative terms. The ansatz for the Lagrangian 
density used by these authors is simply related to our expression for C in Lemma 1. For example, we can add the 
gauge term ^ i^"^^) + ^ (3u^u) + ^ {u2xv) on the right side of (16) and write 

C = v{ut + 6uux + usx) . (33) 

In this context we note that (33) was originally suggested by Atherton and Homsy [ITj and subsequently included 
as an exercise ( 5.37, page 184) in ref. 5. However, it appears that there is no physically founded assumptions in 
writing C in this form except that a Lagrangian may involve its own equation of motion provided one introduces a 
new concept of variational symmetry called the s-equivalence jisl . On the other hand, the treatment presented here is 
based on a fomalism that is specially intended to bring out the reasons why nonloinear evolution equations, as such, 
do not follow from the action principle. 
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The KdV-like equations support bright sohtons while the associated fields have soliton solutions which are dark. 
The coupled set of equations for the bright and dark solitons are Lagrangian although, individually, each of them is 
non-Lagrangian. This observation appears to bring in some similarity with the celebrated work of Bateman Q on the 
dissipative system in particle dynamics. 

The appearance of dark solitons in the solutions of our coupled equations is not a strange phenomenon. In the 
past, dynamics of dark solitons induced by stimulated Raman effect in the optical fiber were explained by the use 
of KdV-Burgers equation jl^. We feel that the inverse variational problem as treated here will be useful to study 
Noether symmetries and even to construct exact solutions of the nonlinear evolution equations [20| . 
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